REGULARITY CRITERIA AND UNIFORM ESTIMATES FOR 
THE BOUSSINESQ SYSTEM WITH THE TEMPERATURE- 
DEPENDENT VISCOSITY AND THERMAL DIFFUSIVITY 



JISHAN FAN, FUCAI LI*, AND GEN NAKAMURA 

Abstract. In this paper we establish some regularity criteria for the 3D 
Boussinesq system with the temperature-dependent viscosity and thermal dif- 
fusivity. We also obtain some uniform estimates for the corresponding 2D case 
when the fluid viscosity coefficient is a positive constant. 



1. Introduction 

In this paper we establish some regularity criteria for the 3D Boussinesq system 
with the temperature-dependent viscosity and thermal diffusivity ( [14,33]): 

dtu + u-Vu + Vn ~ div (^(6')(Vu + Vu'^)) = 6*63, (1.1) 

divw = 0, (1.2) 

dtO + u-Ve- div {k{9)V9) =0, (1.3) 

iu,e)ix,0)^iuo,eo)ix), xeM.^. (1.4) 

Here the unknows u, tt and 9 denote the velocity, pressure and the temperature 
of the fluid, respectively. 63 := (0, 0, 1)"^ is the unit vector in the vertical 0:3- 
direction. The kinematic viscosity ii{0) and the thermal conductivity k{0) are 
generally depending on the temperature 9. Throughout this paper, we assume that 
IJ,{9) and k{9) are smooth functions and satisfy 

< C^^ < fi{9),K{9) <Ci<oo when 16*1 < C2 (1.5) 

for some positive constant Ci > 1 and C2 > 0. Boussinscq type system (1.1)- 
(1.3) are used to model geophysical flows such as atmospheric fronts and ocean 
circulations (see [31,35]). 
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There are many important progress on the Boussinseq system when both k(9) 
and fi{e) are independent oi 9, see [2-4,6-11,15,16,19-22,24,27,28,32] and the 
references cited therein. When k{9) and fi(6) depend on 9, only a few results are 
aviable, see [1,12,13,29,30,37]. In [29,30], Lorca and Boldrini proved that the 
global existence of weak solution with small initial data and the local existence of 
strong solution for general data to the problem (1.1)-(1.4), see also [12, 13]. In [1], 
Abidi obtained the global existence and uniqueness result to the problem (1.1)-(1.3) 
with fc = in some critical spaces when the initial temperature is small. Recently, 
Wang and Zhang [37] obtain the global existence of smooth solutions to the problem 
(1.1)-(1.4) with general initial date in two-dimensional case. However, whether or 
not the smooth solution to the 3D boussinesq system (1.1)- (1.3) with general initial 
data blows up in finite time is still a big open problem. In this paper we provide 
some regularity criteria related to this topic. We will prove 

Theorem 1.1. Let k{9) = 1 and ^i{9) satisfy (1.5) and uo,9q <E H^{M.^) with 
divuo = m M'^. Suppose that there exists a T > such that one of the following 
two conditions holds: 

u e LP(0,T;L«(M3)) with - + - = 1, 3<g<oo, (1.6) 

P q 

Vu e LP{0, T; L9(m3)) with - + - = 2, 3/2 < q < oo. (1.7) 

P q 

Then the solution {u,9) of the problem (1.1)-(1.4) satisfies 

u,9e L°°iO, T; H^{«.^)) f] L^{0, T; H^{M.^)). (1.8) 

Theorem 1.2. Let fi{9) and k{9) satisfy (1.5). LetuQ,9o e i7^(R3) withdivuo = 
in M.'^ . Assume that (1.6) holds with 3 < q < oo, then the solution {u,9) of the 
problem (1.1)-(1.4) satisfies (1.8). 

Remark 1.1. We can obtain the regularity criteria (1.6) or (1.7) for the Cahn- 
Hilliard-Navier- Stokes .system [18, 23]: 

dtu + u- Vu + V-K - V ■ fjL{(l)){Vu + Vu^) = ?yV(/), 
div u = 0, 

dt<j) + u ■ V (j) - Ar/ = 0, 

77 = -A<^ +/'(</>), /(^) ^ i(<^2 _ 1)2^ 

(u, (t>){x, 0) = (uo, 0o)(a;), x g 
Since the proof is similar to Theorem 1.1, we omit the details here. 
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To prove Theorems 1.1 and 1.2, it suffices to obtain a priori estimates under the 
conditions (1.6) or (1.7). Since the equations (1.1) and (1.3) are strongly coupled, 
it is difficult to obtain the estimates of W9 and the Holder continuity of 6, we 
shall use the results in [5,37], Sobolev imbedding, and the ingenious applications 
of Holder and Gagliardo-Nirenberg inequalities to overcome it. 

Next, we consider the Boussinesq system (1.1)-(1.3) in two-dimensional spatial 
space and let the kinematic viscosity fi{6) be a positive constant, i.e, iJ.{6) = e > 
0. In this case the system reads 

Stu + u • Vu + Vtt - eAu = 6*62, (1.9) 

divw = 0, (1.10) 

dt9 + u-Ve-div{K{e)V9)=0, (1.11) 

{u,e){x,0) = iuo,Oo)ix), xeR^. (1.12) 

Here 62 := (0, 1)'^. As before, the thermal conductivity k{6) is assumed to be 
smooth and satisfy 

< Cf 1 < K{e) < Ci < 00 when |6i| < C2 (1.13) 

for some positive constant Ci > 1 and C2 > 0. 

Notice that in [37] Wang and Zhang have established the global existence of 
classical solutions for the problem (1.9)-(1.12) with general initial data and fixed 
e > 0. In this paper we give the uniform-in-e estimates for smooth solutions to the 
problem (1.9)-(1.12). Our result reads 

Theorem 1.3. Let k{9) satisfy (1.13) and e e (0, 1). Let Oq e H^{R^) and uq G 
H^lM."^) with divuo ~ in M.^ . Then for any given T > 0, the solution {u,9) to the 
problem (1.9)-(1.12) satisfies 

l|w||L~(0,T;_f/3(K2)) + |[6'||l~(0,T;H2(R2)) + ||^^||L2(o,T;ff3(K2)) < C, (1-14) 

where C is a positive constant independent of e. 

Finally, considering the following 2D Boussinesq system 

5tu + u- Vu + Vtt- Au = 6'e2, (1.15) 

divu = 0, (1.16) 

546* + u • V6I - ediv (k(6I)V6I) = 0, (1.17) 

(u, e){x, 0) = (uo, 6'o)(a;), x e M^, (1.18) 

we can obtain a similar result on uniform-in-e estimates for smooth solutions to the 
problem (1.15)-(1.18). 
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Theorem 1.4. Let k{9) satisfy (1.13) and e e (0,1). Let uo,^o e H'^{M?) and 
div Mo = m K^. Then for any given T > 0, the solution {u,9) to the problem 
(1.15)-(1.18) satisfies 

||u|iL°°(o,T;ff2(R2)) + ||u||l2(o,T;_h3(k2)) + ||6'||l°°(o,T;H2(r2)) < C, (1-19) 

where C is a positive constant independent of e. 

Remark 1.2. In [7], Chae has proved similar results to Theorems 1.3 and 1.4 when 
k{6) = 1, it seems impossible to apply his method directly to our case. In fact, if we 
follow his arguments we shall encounter the unpleasant terms involved k{9) which 
are out of control. 

We give some comments on the proofs of Theorems 1.3 and 1.4. Since the 
local-in-time well-posedness has been proved in [30], to complete our proofs, we 
only need to prove thea priori estimates (1.14) and (1.19). We shall employ an 
elaborate nonlinear energy method to obtain these desired bounds. More precisely, 
we first use the maximum principle to obtain the L°° estimates of 6. Next, we derive 
an energy estimate based on the L^ energy. Finally, we use Amann's estimates 
on the uniform parabolic equations, the Sobolev imbedding, bilinear commutator 
estimates, logarithmic Sobolev inequality, and Gagliardo-Nirenberg inequalities to 
obtain the desired higher order spatial estimates on 9 and u. 

This paper is organized as follows. In Section 2, we give some notations and 
recall some basic inequalities. Theorems 1.1, 1.2, 1.3, and 1.4 are presented in the 
subsequent four sections. 

2. Preliminary 

In this section, we recall some basic inequalities which shall be used frequently. 

Lemma 2.1 (Logarithmic Sobolev inequality [26]). For all u e i/*(R'') with s > 
1 -|- d/2, there exists a constant C such that the following estimate holds 

\\\7v\\l^ < C{1 + lldivt-llioc + l|curli;||Loo)(l + log(e + \\v\\h^)). (2.1) 

Lemma 2.2 (Gagliardo-Nirenberg inequality [17,34]). Let v e W'''''{R'^)r\L'J{R'^), 
^ l£ l£ oo. Then the following inequalities hold 

\\D'v\\l.<Mo\\D''v\\14v\\1-^, \/0<z<k, (2.2) 

where 

1 i /I k\ . ,1 
p 1\ \r a/ q 
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for all a in the interval 

i 

- < a < 1. 

k 

The constant Mq depending only on d, m, j, q, r and a, with the following exceptional 
case: 

(1) If i = 0, rk < d,q = oo then we make the additional assumption that either 
V tends to zero at infinity or v G L^{M.'^) for some finite q > 0. 

(2) // 1 < r < oo, and k — i — d/r is a non negative integer then (2.2) holds 
only for a satisfying i/k < a < 1. 

We define the operator A := (— A)^/^ via the Fourier transform 
Generally, we define A*/ for s e R as 
For s G K, we define 

11/11^., :=!iAV!u= = ( / levm'd^y^' 

and the homogeneous Sobolev space H'^{W^) := {/ e 5'(M'^) : ||/||^., < oo}. 
Similar, the Sobolev space is equipped with the norm 

Now we recall the following bilinear commutator and the product estimates. 

Lemma 2.3 ( [25]). Let s > 0,1 < p < oo. Assume that f,g e H-'''P{R'^), then 
there exists a constant C , independent of f and g, such that, 

Wifg) - /A-^glU- < v/!U« W-'g\\L^^ + II A^IU-^ II^IU^^), (2.3) 

||A^(/.9)|Up < C{\\f\\L.AWg\\L^. + ||A-7l|L-||g|U,.), (2-4) 
where Pi, p2 € (l,oo] satisfies 

11111 
- = \ = \ . 

p pi qi P2 q2 

Lemma 2.4 ( [36]). Let s > 0,l< p < oo, and f e ij"'P(M'^) n L°°{R'^). Assume 
that F(-) is a smooth function on R with F(0) — 0. Then we have 

||F(/)||^.„ < C{M){1 + ||/|k^)W+i||/||^.., (2.5) 

where the constant C{M) depends on M := sup ||i^'^'^^(t)||^oo . 

fc<M+2,|t|<|l/||i=o 
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In the subsequent sections, we use C (independent of e in Sections 5 and 6) to 
denote the positive constant which may change from hne to hne. We also omit the 
spatial domain R'^ or ^ in the integrals below for simplicity. 

3. Proof of Theorem 1.1 
First, it follows from the maximum principle and (1.3) and (1.4) that 

ll^ilL~(0,T;L~(R3)) < ||0o||l°'(R3) < C. (3.1) 

Multiplying (1.3) by 9, integrating the result over M'^, and using (1.2), we see 
that 

li- I 6'dx 



j e^dx + j \\/e\^dx 0, 



2dt 
which gives 

li^llL°=(0,T;L2(R3)) + ||^||l^(0,T;H1(R3)) < C- (3-2) 

Multiplying (1.1) by u, integrating the result over M? , and using (1.2), (3.1), 
(3.2) and (1.5), we find that 

~ Ju^dx + ^ J \Wu\^dx< J 0e3-udx<\\u\\L2\mL^ <C\\u\\l2, 

which implies 

II"I1l~{0,T;L2(r,3)) + ||u||L2(o,T;ffl(K^)) - 

(I). Assume (1.6) holds. 

Applying the operator di to (1.3), multiplying the result by {di9)^, integrating 
over R^, summing over i, and using (1.2), the Holder and Gagliardo-Nirenberg 
inequalities, we have 
1 d 



4dt 



\Ve\^dx + ^ J \V\Ve\^\^dx <C J |u||V6'nv|V6'nda; 



L 9-2 

|2||l-ai||r7IV7fl|2||l+Q!l 



1 

Hence it holds 

||V0|Uoo(o,T;L4(R3)) < C. (3.4) 

Multiplying (1.1) by — Au, integrating the result over R'^, and using (1.2), (3.4), 
the Holder and Gagliardo-Nirenberg inequalities, we derive that 
1 d /■ ,„ ,9 . 1 



2 dt 



J \\7u\^dx + ^J \Au\^dx< J {u ■ \/)u ■ Audx + C J |V6'||Vw||Au|da; 



<||u||l<!1|Vu|| j2^\\Au\\l2 + C||V6i!|i4j|Vui|L4||Aui|L2 

L 9-2 
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<C!|u||L,||Vu||^n|Au||^+"^ + C||Vu!U4||Au||i2 
<2^l|A^||i.+C{h||^,||V^||i. + ||Vu||i.}, 
which gives 

lkllL°=(0,T;Hi(R3)) + l|w||L2(0,T;ff2(R3)) < C. (3.5) 

Applying the operator A to (1.3), multiplying the result by A0, integrating over 
M'^, and using (1.2), (3.5), the Holder and Gagliardo-Nirenberg inequalities, we 
obtain 

j \M\^dx + j \VM\'^dx^^ j d^{u-ye) -d.Aedx 

i 

<{\\uhs\\Ae\\Ls + ||V7i|U.||v0|U^.)||vA0|U2 

<C{\\A9\\LS + \\W9\\L^)\\WAd\\L. 



<C\\A0\\Y.'\\VAe\\%' < -\\WAe\\l.+C\\A0\\l,, 



1 

2 

which leads to 



II^IIl~(0,T;H2(r3)) + ||^||L2(o,T;ff3(K3)) < C. (3.6) 

Applying the operator A to (1.1), multiplying the result by Au, integrating over 
^ and using (1.2), (1.5), (3.1), (2.5), and (3.5), we arrive at 

~ J \Au\'dx+^ J \VAu\'dx 
<YlJ ^^{u■Vu) ■ d^Audx + C j {\A^i{e)\\Vu\ + \V^i{9)\\V'^u\)\VAu\dx 

i 

<(||ii||i6||Au||i3 + \\\/u\\l2\\Vu\\l'^)\\VAu\\l2 

+ C{\\Ae\\L3\\Vu\\LS + \\Ve\\L^\\Au\\L2)\\VAu\\L2 
<C{\\Au\\ls + \\Vu\\l^)\\VAu\\l2 

+ Ci\\Ae\\L3 + \\Ve\\L^.)\\Au\\L2\\VAu\\L2 



<-\\WAu\\l2 + C\\Au\\l2 + C{\\Ae\\l, + \\W0\\U)\\Au\\l2, 



which implies 

lkllL°=(0,T;H2(R3)) + 11^11^2(0, T;ff3(R3)) < C. (3.7) 

Thus (1.8) holds. 

(II). Assume (1.7) holds. 

Applying the operator di to (1.3), multiplying the result by {di6)^, integrating 
over M"^, summing over i, and using (1.2), the Holder and Gagliardo-Nirenberg 
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inequalities, we obtain 

1 d /",„„, 4 , . 3 
Idi 



<c J |Vw||v6i|2 . \\/e\^dx 



a2 



<c\\vu\\L4\'^d\Y ^ 

<^\\'^mYL^ + c\\vu\\u\\yo\YL^ (p = Y 

which proves (3.4). 

Multiplying (1.1) by —Aw, integrating the result over M.^, and using (1.2), (3.4), 
the Holder and Gagliardo-Nirenberg inequalities, we get 

<YlJ u^^^u^]udx + C J | V6i| | Vu| | Au|(ia; 

<-^ / dju.d^udjudx + C\\\/e\\L4\\\/u\\L4Au\\L2 
i 

<C||Vu||i,||VM||2^ +C\\\/u\\l4Au\\l2 
<C\\Vu\\L4Vu\\l^t"'^\\Au\\lT + C\\Vu\\l4Au\\l2 
<^||A^.||i. + C{\\Vu\\l4Vu\\l2 + llVHIi.}, 
which implies (3.5). 

Noticing that the calculations for (3.6) and (3.7) still hold in this case, we hence 
complete the proof of Theorem 1.1. 



4. Proof of Theorem 1.2 
First, by the maximum principle, it is easy to prove that 

I1^!1l='(o,T;L°°(r3)) < C- (4.1) 

Noticing the condition (1.5), we still have (3.2) and (3.3). 
Similar to the Proposition 4.1 in [37], we can prove 

Proposition 4.1. Let u satisfy (1-6) and divu = in R"^ x (0,oo). Assume that 
e e L°°{0,T;L'^) n L'^{0,T;H^) is a weak solution to (1.3) and (1.4). Then there 
exists aae {0, 1) such that 9 e C"([0,r] x R^). 
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Proof. Most of the calculations are as same as that in [37], the only difference is 
the calculations of the following term 

rt2 



a2 



u ■ V77 rjdxdt 



t2 



< / \\u\\L4vOk\\ ^||0feV?y||i2di 
Jt, Li-'-' 



ti 

t-2 

2 \\„a l|2(l-"i)||Y7/'„fl M|2ai 7, , r<\\a Y7„l|2 



1 — ai 
-=(ti,t2;L2(R3)) 

+ C'e||VMfe)||i2(t^_4^.i2(R3)) + C'||6lfeV77|1^2(t^_t^.i2(H3)) 

for any < e < 1. □ 
Now, using an estimate of the gradient of the solution to the parabolic equation 
dtO - div {K{e)Vd) = -div {uO), 

we have (see [5]) 

1|V^1Ilp(0,T;L9(R3)) <C'lk^llLP(0,T;L9(R3)) + C 

<C||6'||ioo(o,T;L~)ll"llLP(0,T;L'!(R3)) + C* 

<C'II^IIlp(o,T;L9(r3)) + C"- (4.2) 

Multiplying (1.1) by — Au, integrating the result over M? , and using (1.2), (4.2), 
the Holder and Gagliardo-Nirenberg inequalities, we have 

i| j \Vu\^dx + ^ J \Au\'dx 

< j {u ■ S/)u ■ Audx + C j |V6l||Vu||AM|da; 

<||u||L9i|Vu|| ^||Ait||i2 + C||V6'||L,||Vit|| ^||Ait||i2 
<C{\\u\\li + ||V0||LOI|V«||i2"M|Aw||it"^ 
<^||A^i||i2+C(||«||i, + ||V^?||i,)||V«|li2, 
which gives (3.5). 

Applying the operator A to (1.3), multiplying the result by A^, integrating 
over R^, and using (1.2), (3.5), (2.5), (4.2), the Holder and Gagliardo-Nirenberg 
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inequalities, we have 
1 d 



2 dt 



j \M\^dx + ^ j |VA6ipdx 



<(||i/,|U6||A0|U3 + ||V^.|U2||V0||l=.)||VA0|| +q|V0|U,||A(?||^^"i||VA0||i+"i 

<c(||A0|U3 + i|v^^iU^)||vA0|U. + c\\ve\\u\M\\l.. + ^llvA^iii^ 
<^\\vM\\i. + c{\\M\\h + \\ve\\UM\\i.], 

which imphes (3.6). 

Noticing the calculations for (3.7) still hold here, we hence complete the proof 
of Theorem 1.2. 



5. Proof of Theorem 1.3 

To complete our proof of Theorem 1.3, we only need to prove a priori estimates 
(1.14). By the maximum principle, it follows from (1.11) and (1.12) that 

II^I|l°°(o,T;L°=(r2)) < ||0o||l~(r2) < C. (5.1) 

Multiplying (1.11) by 9, integrating the result over R^, and using (1.10), (1.13), 
and (5.1), we see that 

1 d 
2di 
which gives 

I1^IIl = (0,T;L2(r2)) + ||0||l2(O,T;H1(R2)) < C. (5.2) 

Multiplying (1.9) by m, integrating the result over R^, and using (1.10) and (5.2), 
wc find that 



'^dx+ / K{e)\V0\'^dx = 0, 



J u^dx + t j \Vu\'^dx = j 9e2udx 



< \\9\\L4^i\\L2<C\\u\\L2, 



which yields 

I|w|1l°=(0,T;L2(r2)) + \/el|M||L2(o^T;Hi(K^)) - (^■^) 

Multiplying (1.9) by —Au, integrating the result over M^, and using (1.10), (5.2) 
and the fact that 

J{u ■ V)ii • Audx ^ 0, 

we get 

ld_ 
2dt 



\Vu\^dx + e J \Au\^dx< J |V6'||Vu|dx < |lV6i|li2|lVu|U2, 
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which leads to 

II"I|l°°(0,T;H1(K2)) + \/e||^illL2(o,T;_H-2(R2)) < C. (5.4) 

Then by the same calculations as that in [37] , we have 

||^||c°(E3x[0,T]) < C (5.5) 

for some a G (0, 1) independent of e > 0. 

Now, using an estimate of the gradient of the solution of the parabolic equation 

dtO - div {n{e)Ve) = -div {uO), 

we obtain (see [5]) 

I|V6'|1l4(o .r.^4(R2)) <C + ||u6'j|i4(o^T;L-i(]R2)) 

<C + ||u||l4(o,T;L1(r2))||^!1l°=(o,t;L=(k2)) < C. (5.6) 

Applying the operator A to (1.9), multiplying the result by A0, integrating over 
M^, and using (1.10), (5.4), the Holder and Gagliardo-Nirenberg inequalities, we 
derive 



< h^i^ ■ ^0)d,Mdx + C||V0|U4||A0|U4|1VA0|U2 



<c{\\u\\LA\m\L- + \\Vu\\LA\'^e\\L^)\\vAe\\L- + c\\ve\\LA\mLA\'^^o\\L^. 

<C{\\Ae\\L^ + ||V0||l~)||VA0|U2 + C||V0|U4||A0|U4||VA0|li2 
<C(||A(?||^/'||VA0||^/' + ||V0||'/'||VA0||^//)||VA^^||i2 

+ C||V0|U4||A0||^/'||VA(?||^/' 
<^||VA0||i. + C{\m\l._ + \\V9\\l, + |lV0|li.|lA0|li.}, 
which implies 

II^'I|l°=(0,T;H2(r2-)) + ||^||L2(o,T;if3(K2)) < C. (5.7) 

Here we have used the estimate (see (2.5)): 

||A«;(6')||lp < CIIAe'lliP with l<p<cx). (5.8) 

Let uj = curlu = diU2 — d^ux for u = {ui,U2)^ ■ Taking curl to (1.9) and using 
(1.10), we infer that 

atw + u- Vw-eAw^^ai^. (5.9) 
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Multiplying (5.9) by integrating the result over M^, and using (1.10) and 

(5.7), we have 

1 d 



q dt 

which gives 



\^Vdx < / die\uj\''-^ujdx < ||9i6ij|L9||wf 



j^Ml.<C\\V0\\l., 



and thus 

II^IIl. < IIc^oHl. + c / \\ve\\L.dx. (5.10) 

Taking q — > +00 in (5.10), we obtain 

II^^I|l°=(o,T;L~(r2)) < C. (5.11) 

Applying the operator A"^ to (1.9), multiplying the result by A'^m, integrating 
over R2, and using (1.10), (5.7), (5.11), (2.3), and (2.1), we conclude that 

~ J \A^u\^dx <- J {A^u ■ Vu) - uVA^u)A^udx + J A^Oe^) ■ A^udx 

<C\\Vu\\lo.\\A\i\\12 + C\\A^e\\L2\\A^u\\L2 

<C(1 + log(e + \\A\Um\A\\\l2 + C\\AmL4^\U2, 

which leads to 

II^IIl=-(o,T;H3(ji;2)) < C. (5.12) 
This completes the proof of Theorem 1.3. 

6. Proof of Theorem 1.4 

To complete our proof of Theorem 1.4, we only need to prove a priori estimates 
(1.19). First, the estimate (5.1) still holds. Similarly to (5.2) and (5.3), we have 

ll^'ilL°°(0,T;L2(R2)) + \/e||^||L2(o,T;Hi(R2)) < C\ (6.1) 

||'«||l~(0,T;L2(R2)) + ||u||i2(0,T;_ffl(R2)) < C. (6.2) 

Multiplying (1.15) by —Am, integrating the result over R^, and using (1.16), 
(6.1) and the fact that 

{u ■ V)u • Audx = 0, 



\Vu\^dx + / |Au|^dx = - / 9e2Audx < \\e\\L2\\Au\\L2 < C\\Au\\ 



we find that 
ld_ 

which implies 

I|w||l=^(0,T;_H1(K^)) + ll^l|L2(0,T;ff2(R2)) < C. (6.3) 
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By the Sobolev embedding theorem, it is easy to find that 

dtu- Au + Vtt = f, f ■.^ee2-u-Vu<^ L^{Q,T;LP{M?)), (6.4) 

for any p € (2,oo). Therefore, it fohows from (6.4) and the regularity theory of 
Stokes equation that 

II"I|l2(o,T;W2,p(h2)) < C|1/|1l2(o,t;Lp(k2)) + C <C. (6.5) 

And thus 

I|Vu||l2(0,T;L~(R2)) < C'|Kt||L2(0,T;W2.P(R2)) < C. (6.6) 

Let ^ 

K{e) := / ^(Ode 



We deduce from (1.17) that 

dtK{e) + u ■ Vn{e) - eK{0)ATi{e) = 0. (6.7) 

Muhiplying (6.7) by ~Ak{9), integrating the resuh over R^, and using (1.16) and 
(6.6), we derive 

14- I |VK(0)|2dT + e / K{e)\AK{e)\''dx 



2dt 



2 

L2, 



= j u\/K{e)AK{e)dx = - j J2djuAK{0)djK{e)dx < ||Vu|1loo|ivk( 

which imphes 

ll^^l|L=^(o,T;_f/i(R2)) + V^\\Ak{9)\\ ^Qrp.j^2(^^2)^ < C. (6.8) 
Noting that (see (2.5)) 

Vel|A6'||L2(o,T;L2(R2)) =Ve|lA(K~-^ ° '5(^'))||l2(o,t;L2(r2)) 

<V^C\\AK{9)\\mo,T;L^(«.^)) <C. (6.9) 
By the Gaghardo-Nirenberg inequahty, it foUows from (6.8) and (6.9) that 

e||ve|li4(o.T;L4(R2)) =e r \\ve\\\,dt<ce r wvewumil.dt 

JQ Jo 

<Ce [ \\Ae\\l2dt<C. (6.10) 

JQ 

Applying the operator di to (1.17) gives 

dtd.e + uvd.e + ^^u ■ ve - ediv {K{e)vdie + k'{e)d,eve} = o. (6.11) 

Multiplying (6.11) by (diO)^, integrating the result over K.^, summing over i, and 
using (1.16), (6.6) and (6.10), we obtain 
1 d f ,1 . e 



4dt 
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<Ce J \ve\-\\7e\^\\/\\/e\^\dx + c\\Vu\\L^ J \ye\*dx 

<Ce||V^?iU4||(V0)2||^4i|V|V^^|2||i2 + C\\Vu\\l^ J 
<Ce\\\/9\\L4\\/9\^\\'/^ -WVlVemli^ + C\\\7u\\l^ J \V6\Ux 

<^\Myo\YL-^ + Ce\\ve\\U\\yo\'\\h + c\\VuU^ J \ve\Ux, 

which gives 

||V0|Uoo(o,T;L4(R2)) < C. (6.12) 

Applying the operator A to (1.17), multiplying the result by A9, integrating 
over M^, and using (2.5), (6.12), (5.8), (6.5), and (6.6), we conclude that 



1 d 
2di 



J \A9\^dx + -^ J \VA9\^dx 

<- J {A{u ■ \79) - uVA9)A9dx + Ce||V6i||i4|lAfc|li4|lVA6'||i2 
<C(||A^/,lU4||V^?|U4 + ||Vu|U^.||A0|U2)||A0|U. +Ce||A0|U4||VA0!U2 
<G{\\AuU. + ||Vu|U~|lA0|U.)||A0|U. + ^llVA^IIi, + Ce\\A9\\l., 

which leads to 

II^!Il°=(o,T;H2(r2)) < C. (6.13) 

Finally, applying the operator A to (1.15), multiplying the result by Au, inte- 
grating over M^, and using (1.16), (6.13), and (6.6), we arrive at 

j \Au\^dx + j \VAu\^dx 

A9e.2 ■ Audx - A{u ■ Vu) • Audx 



<\\A9U2\\Au\\l2+C\\Vu\\l^\\Au\\1. 
<C\\Au\\l2+C\\^u\\l^\\Au\\1,, 



which gives 



l|w||L~(0,r;ff2(R2)) + ||u||i2(0,T:fl'3(R2)) < C. 

This completes the proof of Theorem 1.4. 
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